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We calculate the ground state energy and phase diagram of the quasi-two-dimensional (quasi-2D) dipolar 
Fermi gas by variationally mapping the system to its strictly 2D counterpart for which accurate quantum Monte- 
Carlo (QMC) results are available. We further demonstrate the accuracy and flexibility of our variational tech- 
nique by mapping the 2D dipolar Fermi gas to the 2D electron gas and showing that the resulting variational 
estimates closely match the QMC results. Our analysis suggests that the strongly-correlated Fermi liquid wave- 
functions of different 2D Fermi systems share salient universal features. We show that inclusion of correlations 
significantly modify the dependence of the properties of the quasi-2D dipolar Fermi gas on the thickness of the 
layer compared to mean-field predictions, and calculate the frequency of its monopole oscillations in isotropic 
traps as an experimental demonstration of correlation effects. 



Dipolar quantum gases have been the subject of much ex- 
perimental and theoretical interest in the recent years. The 
long-range and anisotropic nature of dipole-dipole interac- 
tions can be exploited to realize a wide range of genuinely 
new quantum many -body phases | IJ. Experiments with dipo- 
lar bosons have made rapid progress |I2}S1> with fermionic 
polar molecules and rare earth atoms emerging as particularly 
promising candidates for realizing strongly interacting dipo- 
lar Fermi gases owing to their large permanent dipole mo- 
ments |5-9|. The collisional stability of trapped dipolar gases 
requires suppression of the attractive dipole-dipole (head-to- 
tail) scattering. This is achieved by tightly confining the gas 
in a thin "quasi-two-dimensional" layer and polarizing the 
dipoles perpendicular to it. The resulting effective dipole- 
dipole interaction is purely repulsive, has a long-range r^^ 
tail, and produces a Fermi liquid (FL) state with novel proper- 
ties LIOJ . 

Although the theoretical investigation of dipolar fermions 
began more than a decade ago, the vast majority of obtained 
results rely on mean-field theory 1 1 1 1. The natural next step is 
to investigate the role of correlations. Recently, Matveeva and 
Giorgini have caiTied out a fixed-node diffusion Monte-Carlo 
(FN-DMC) analysis of strictly 2D dipolar Fermi gas (2DDFG) 
at zero temperature and have shown a direct FL to triangu- 
lar Wigner crystal (WC) transition in the strongly interacting 
regime. The mean-field theory predicts a density-wave insta- 
bility as well lfT3l . albeit at an interaction strength almost 20 
times smaller than the QMC result. 

In experiments with 2D dipolar fermions, the transverse 
confinement potential has a finite strength and the quasi-2D 
layer has a finite thickness as a result. The inclusion of this 
effect in realistic calculations is of particular importance here 
due to the strong anisotropy of dipole-dipole interactions. In 
this Letter, we use a simple variational method introduced 
earlier by two of the authors fT4',T31 to incorporate the effect 
of finite transverse confinement in the ground state energy 
and the phase diagram of the quasi-2D dipolar Fermi gas 
(qDFG). This is achieved by using the variational principle 




FIG. 1. (Color online) (a) The energy of quasi-2D dipolar Fermi 
gas (qDFG) as a function of dipole-dipole interaction strength gd = 
kpCLd (defined in the main text) and the layer width in the 
units of £o = h^Trn/m. The upper surface and the accompanying 
black dashed lines show the result from Hartree-Fock theory while 
the lower surface and the black solid lines represent the variation- 
ally obtained result using 2DDFG ground state wavefunctions. The 
solid red lines show the variationally obtained energies using 2DEG 
ground state wavefunctions. The weak -dependence of the vari- 
ationally obtained results is due to correlations, (b) The variation- 
ally obtained phase diagram of qDFG. The MS region indicates the 
breakdown of the single-subband approximation (cf. main text for 
details). The dashed lines indicate the lower and upper uncertainty 
bounds for the WC transition line. 



to approximately map the ground states of qDFG to a known 
system such as the two-dimensional electron gas (2DEG) 
or 2DDFG. In the first stage of the calculations, we assess 
the reliability of the variational technique by estimating the 
ground state energy and the WC transition point of 2DDFG 
using the ground states of the two-dimensional electron gas 
(2DEG) as tiial wavefunctions. We show that the obtained 
variational estimates are remarkably close to the QMC 
results, capturing more than 95% of the coiTelation energy in 
the strongly interacting regime (see Fig. |2^) and predicting 
the coiTect location of the WC transition within the accuracy 
of QMC calculations. Although 2DEG and 2DDFG are 
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described by very different microscopic interactions and 
undergo crystallization in different regimes (low and high 
densities respectively), this analysis implies that these two 
seemingly different systems exhibit a very similar behavior 
in the strongly correlated regime even from a quantitative 
vantage point. In the second stage of calculations, we 
obtain the ground stage energy and phase diagram of qDFG 
based on 2DDFG wavefunctions (see Fig. [T]i. We show 
that while mean-field theory predicts a strong dependence 
of the properties of qDFG on the transverse confinement 
width ttz, in actuality much of the short -ranged details of the 
effective interaction are masked by the exchange-correlation 
hole, resulting in a much weaker -dependence. We finally 
discuss the experimental signatures of strong correlations in 
experiments with ultracold dipolar gases. We calculate the 
frequency of collective monopole (breathing) oscillations of 
qDFG in harmonic traps and show that correlations produce 
an observable effect on them. 

The model. - We consider a gas of spinless dipolar fermions 
at zero temperature in an oblate and axially-symmetric confin- 
ing potential J7trap(r) = mcjQ(a;^ + J/^)/2 + ma;^z^/2, such 
that Wz ^ wq. We focus on the quasi-2D limit hiOz ^ ep (^f 
is the Fermi energy) where only the lowest transverse subband 
is occupied. In this limit, the system can be described as a 2D 
Fermi gas with the following two-body interaction potential: 

V2fW= /d^dz'|0oWn<^o(^')Pv|Pp(r,z-z'), (D 

where V^°(r, z - z') = Z?^ (|r|2 _ 3^2-) /|j.|5 jj^g dipole- 
dopole interaction in 3D space (D is the dipole moment 
and we have assumed that the dipoles are aligned along 
the z-axis) and 4>oiz) = ^^^^^^ /{y^ az)^ is the trans- 
verse wavefunction of particles in the lowest subband (a^ = 
[h/{mujz)]~^^^ is the transverse oscillator length and m is 
the mass of a single particle). The analytical expression for 
V^p°(r) is given in the Supplemental Materials (SM) §1. We 

find that vSp°(r) w [2D^ /{V2^al)]\n{az/r) + 0(1) for 

r < a, and V5p°(r) « D^r^ - 9a2/(2r5) + 0{r-^) 
for large r. The single-subband limit is valid as long as 
/i < hujz, where fj, is the chemical potential and hujz is the 
non-interacting band gap fT9l. 

In the homogeneous case (no in-plane trap potential), the 
ground state can be labeled by two dimensionless parameters, 
gd = kpGd, and = \fnaz- Here, kp = \J ^-kji is the 
Fermi momentum, n is the 2D density and ad = mD^ jf? 
is the "dipolar length". Here, we confine our attention to the 
regime ^ 1, where the single-subband approximation 
is applicable. In the presence of the trap potential, we treat 
the problem in the local density approximation (LDA) and 
parametrize the system with gd = kp^^ad and = ^/n^az, 
where uq = {2N)i muio / {2nh) is the density at the center of 
the trap for a non-interacting Fermi gas, and fc^ o = V^^^rn^. 

The variational mapping method. - Variational methods 



are known to work well if the class of trial states is suffi- 
ciently flexible. Having decent variational wavefunctions at 
one's disposal can obviate time-consuming numerical calcu- 
lations for many practical purposes. The states we consider 
here are the ground states of a repulsively interacting Fermi 
system, such as 2DEG or 2DDFG, whose ground states are 
accurately known. Depending on the strength of interactions, 
these states span a wide gamut, from a free Fermi gas to a 
strongly-correlated, crystal-like liquid. This choice guaran- 
tees that our variational energy estimate is asymptotically ex- 
act at both weak and strong coupling. We may therefore ex- 
pect it to be rather accurate everywhere in between as well. 
In fact, in the intermediate strongly-correlated regime, the 
short-range details of interactions are masked by the local- 
ized exchange-correlation hole while the long-range tail is 
screened, making the wavefunctions of different systems bear 
strong resemblance to each other The practical importance of 
this technique is the great simplicity and speed that it offers: 
by "recycling" what we already know about one strongly cor- 
related system, we come up with decent variational estimates 
for another system at a small fraction of the computation time 
and effort that is usually required. 

The procedure is formally described as follows: let l^*^)"^, 
SkI^*!^'*] g^ifi^^) be the normalized ground-state 
wavefunction, kinetic energy density and the pair-distribution 
function of a given system "A" at a fixed density n. Here, 
is a dimensionless quantity that parameterizes the strength 
of interactions in "A" (e.g. in case of 2DEG, ^ can be taken 
as the Wigner-Seitz radius r^). We propose {I'l'^)'^} as a 
family of variational wavefunctions for a second system "B". 
Let e^I'^'^var] and ej^j[n; ^var] be the kinetic and interac- 
tion energy densities of "B" obtained by choosing 
as a trial wavefunction. Since the kinetic energy opera- 
tor tC = — X^i ^^^f /(2"^) is identical for both systems, 
£g[n;evar] = ^(*«™J/C|4'j_)^ = e^[n-U,]. The inter- 
action energy density of "B" in the same trial state can be 
calculated using the "A" pair-distribution function: 

efnt [n; Cvar] = J ^ 3^ (r; Cvar) v"" (r; C^) 2^r dr, (2) 

where v^{r; is the two-body interaction potential of "B" 
and f ^ is a dimensionless measure of interaction strengths in 
"B". Here we have also assumed, for concreteness, that both 
systems are spin-polarized 2D fermions with isotropic inter- 
actions. Minimizing e-^[n;^var] = ei[?i;Cvai] + ej^^tin; Cvar] 
with respect to ^var, we obtain (1) a variational upper bound 
for the ground state energy of "B", and (2) a mapping 
(/) : ^■^ — > that associates the ground states of "A" 
with the (approximate) ground states of "B". The mapping 
also provides a heuristic method for obtaining the phase 
transitions of "B": if a phase transition occurs at for "A", 
the mapping function predicts the same transition in "B" 
at provided that the inverse mapping exists. The 

accuracy of this simple scheme for predicting phase transition 
lies at the heart of our central premise, i.e. the universality of 
strongly-correlated wavefunctions near the FL-WC transition. 
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We refer to this scheme as variational mapping, and the 
results obtained by mapping to "A" as VMa. Note that only 
the knowledge of 5"*(r;^^) and ^xl'^iC^] required here, 
and not the complete wavefunction. Furthermore, if the 
microscopic interaction is a power-law function, the kinetic 
energy density can be extracted from the total energy using 
the Virial theorem. 

Variational mapping of 2DDFG to 2DEG.- In order to 
assess the accuracy of this method, we calculate the energy of 
the 2DDFG using the ground states of the 2DEG and compare 
to the FN-DMC results. Since our model here is composed of 
a single hyperfine state, we use ferromagnetic 2DEG states. 
We use the energies reported in Ref. |21 1 and the analytical 
representation of the 2DEG pair-distribution function pro- 
vided in Ref. Il22l . Fig. |2^ shows the variationally obtained 
energies (green dashed lines) along with the FN-DMC result 
from Ref. 1 12 1 (solid lines). We find that the VM2DEG scheme 
remarkably captures more than 95% of the correlation energy 
of 2DDFG in the strongly interacting regime (g^ > 20). The 
transition from the ferromagnetic liquid to the WC phase 
occurs at = 28 ± 3 for the 2DEG |23 | (r^ is the Wigner- 
Seitz radius), which variationally maps to 5™^c ~ 28 ± 4. 
This is again remarkably close to the FN-DMC prediction 
5d,wc = 25 ± 3 lfT2l . One can notice that the accuracy of 
the variational method improves for stronger interactions 
(see the inset plot of Fig. |2^). We associate this with the 
growth and localization of the exchange-correlation hole 
in the strongly correlated regime [241, which masks the 
short-range details of interactions better and also effectively 
screens its long-range tail. The difference between the exact 
and variationally obtained pair-distribution functions, shown 
in Fig.|2|5 (rightmost plots), is barely noticeable for large gd- 

Variational mapping ofqDFG to 2DDFG. - Here, we take 
a much smaller step compared to the previous section and cal- 
culate the energy of the qDFG using the DFG ground states. 
We naturally expect the accuracy of the variational mapping 
procedure to be even better in this case, given that the 2D and 
quasi-2D effective dipole-dipole interactions have identical 
long-wavelength limits and only differ in their short-range de- 
tails. We use the 2DDFG energies and pair-distribution func- 
tions reported in Ref. [il2|. The details of the calculations are 
given in SM §2. Here, we also calculate the energy of qDFG 
using 2DEG ground states as a control variational estimate. 
Fig.[T^ shows the obtained energy as a function of gd and 
along with the Hartree-Fock result. The FL-WC phase dia- 
gram obtained from the mapping function is shown in Fig.[l]5. 
We notice that VM2DEG (red lines) and VMdfg (black lines, 
lower surface) energies are very close and follow the same 
trend consistently. We also find that VM2DEG -predicted en- 
ergies are consistently slightly larger than VMdfg energies, 
which simply indicates that DFG wavefunctions constitute a 
slightly better family of trial wavefunctions for qDFG. An in- 
teresting consequence of correlations is the essentially differ- 
ent dependence they produce for the ground state energy on 
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FIG. 2. (Color online) (a) The correlation (Ec) kinetic (ekin), and total 
(etot) energy densities of DFG in the units ehf = (ft'^7rn/m)[l + 
1283d/(457r)] (cf. SM §1 or Ref. (Til)- The red triangles are the 
QMC results from Ref. f 121 and the black solid lines are the analytic 
fits (provided in SM §2 for reference). The green dashed lines are 
our variational results based on 2DEG wavefunctions (VM2deg)- The 
cyan and gray vertical bars represent WC transition as predicted by 
QMC (ga = 25 ± 3) and VM2DEG (gd = 29 ± 4) and coincide within 
the accuracy of the calculations. The inset plot shows the ratio of 
the correlation energy captured by VMideg and QMC as a function 
of gd- (b) The pair-distribution function g{r) of DFG obtained using 
QMC (solid lines) and VM2DEG (dashed lines). The black, green, 
blue and red lines correspond \.o gd = Q (Hartree-Fock), 2, 10, 20. 

a^, as compared to mean-field theory. Indeed, Hartree-Fock 
theory predicts a linear dependence of the total energy on 
for small (see the black dashed lines in Fig. [T^; the ex- 
act Hartree-Fock energy expression is also given in SM §1), 
while VMdfg (and VM2deg) strongly suggest a quadratic Oz- 
dependence (see the solid black/red lines in Fig. [T^). In fact, 
we find that the VMdfg energies can be parametrized to an 
excellent approximation as: 

£qDFG(5'd, a^) « EDFclffd) (co + cigd) al, (3) 

m 

in the parameter regime 5 < < 30 and \/naz < 0.1. Here, 
Co = 12.8 and ci = 1.45 are obtained by fitting to the VMdfg 
results, and £DFG(5d) is given in Ref. lfT2]| and SM §2. The 
quadratic dependence of eqDFG on Gz can be understood in 
simple terms: for large gd, the short-range part of the effective 
dipole-dipole interaction is masked by the correlation hole 
and the a^-dependence of the energy results from the leading 
correction to the interaction in the large-r limit, which is 
^ al/r^ [cf. the discussion after Eq. The FL-WC 

transition line can be parametrized as gd,wc ~ 25 + C2na1„ 
where C2 — 4.82 x 10^. In contrast, the Hartree-Fock theory 
spuriously predicts a linear -dependence for the FL stability 
phase boundary ifTSl . 

Experimental observation of correlation effects. - Among 
the important advantages of experiments with ultracold atoms 
and molecules as compared to traditional condensed matter 
systems are their cleanness and their significantly slower dy- 
namics. These features allow direct measurement of the equa- 
tion of state 1 16 1 and the collective modes with a remarkably 
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mass and momentum currents, i.e. 








FIG. 3. (Color online) The monopole oscillation frequency flmon 
of qDFG in an isotropic trap at T — 0. The yellow surface and 
red lines show the mean-field (Hartree-Fock) result. The blue sur- 
face and lines, and the black lines show the lower and upper bounds 
to Qmon once correlation are included using the VMdfg method. 
The bounds are obtained from hydrodynamic theory (lower) and sum 
rules. Inclusion of correlations changes the -dependence of S7mon 
significantly. 



high accuracy ifTTl . We showed that correlations have a pro- 
found effect on the properties of qDFG in the previous sec- 
tion. Here, we show that these effects can be directly ob- 
served by measuring collective oscillations of qDFG in op- 
tical or magnetic traps. Precise measurements of such have 
been successfully utilized to study the physics of the BEC- 
BCS crossover ifTTl . Here, we study qDFG in isotropic traps 
and calculate the frequency of monopole (breathing) oscilla- 
tions, Onion, about the equilibrium state. 

Trapped qDFG. - The LDA equilibrium state of the trapped 
gas can be obtained from the balance of the trap restoring 
force and the internal pressure of the gas, i.e. drP[ncq{r)] + 
neq(r) 9rJ7trap(r) — 0. Here, noq(r) is the local equi- 
librium density and P[n] is the LDA pressure functional. 
At zero temperature, the pressure is given by P[n] = 
dn' n' d{ijL[n'] ) / 9„/ , where fi[n] — d{ne) jdn is the chem- 
ical potential. The equilibiium condition yields: 



dn, 



cq 



= -r7iW(,K[rieqjn^qr, 



(4) 



where = n^^[9^(n£[n])/(9^n]^^ is the compressibility, 
which can be calculated from the energies presented above. 
The equilibrium density profile is obtained by solving Eq. Q 
under the global particle number constraint, J d^rncq(r) — 
N. Correlations reduce the pressure at a given density, mak- 
ing the gas more compressible and resulting in a consistently 
smaller equilibrium radius of the trapped gas as compared to 
mean-field theory. 

An exact treatment of collective oscillations at T = re- 
quires solving the Landau kinetic equation, which is not fea- 
sible without the knowledge of the Landau parameters. Nev- 
ertheless, a lower bound to the frequency of collective oscil- 
lations can be found by resorting to the ideal hydrodynamic 
approximation. At zero temperature, the hydrodynamical de- 
scription of the gas is given by the conservation laws for the 



m{dt'v 



dtn 



dr{nv) = and 
where v is the 



macroscopic velocity field. Linearizing these equations about 
the equilibrium state and solving for Sn = n — n^q gives 



dt dn - 



(5) 



The absence of dynamical Fermi surface deformations in 
the hydrodynamic approximation (which strictly increase the 
sound velocity in the case of isotropic repulsive interac- 
tions), implies that the obtained collective excitation ener- 
gies are strictly lower than the exact values 12514271 . On the 
other hand, a stiict upper bound is given by (Srimon)^ < 



''M 



jm 



(1) 



where mVj and m'\l are the first and third mo- 



,(3) 



ments of the monopole response function and can be evalu- 
ated using sum rules 1 28 1 . The numerical method for solving 
Eq. (|5]) and calculating the required moments are discussed in 
detail in SM §3. The mean-field analysis of collective modes 
in traps has appeared in several previous works ifTTHTSl . Here, 
we obtain 51i„on by solving the Boltzmann-Vlasov equation 
using the numerically exact method given in Ref. fTS]- 

Fig.|3]shows the obtained results. We note that the different 
a^-dependence of the mean-field vs. correlated theory also 
persists in rimon- For small values of a^, correlations have 
a tendency to decrease the frequency of oscillations as com- 
pared to the prediction of the mean-field theory. This scenario 
is reversed, however, as is increased. One can understand 
this reversal as a consequence of the weak dependence 
of the correlated theory on the width of the layer and the 
erroneously strong dependence of the mean-field theory. We 
note that while strictly speaking the results presented here 
correspond to the T = limit, the correlation effects are 
expected to persist as long as T < Tp. In particular, the 
quadratic dependence of 57nion on the width of the layer only 
relies on the presence of the correlation hole and is expected 
to persist up to T ^ Tp. Suitably large dipolar interactions 
for observing such strong correlation effects are expected to 
be achievable in experiments with fermionic polar molecules 
such as KRb S and NaK E). 

Summary and outlook. - We have explored the idea of 
variational mapping between the ground states of different 
Fermi systems as a simple method for obtaining ground state 
energies of strongly correlated fermionic systems. We have 
found strong quantitative evidence for the accuracy of this 
simple variational scheme, which suggests a universality 
in strongly-correlated Fermi liquid wavefunctions. We 
have calculated the energy and the phase diagram of qDFG 
using this method and shown that coiTelations result in very 
significant connections to the mean-field picture. Finally, 
we have calculated the frequency of monopole oscillations 
of qDFG as a simple experimental proposal for observing 
the correlation effects. As a final remark, we note that the 
variational mapping method may be extended to superfluid 
and multi-component Fermi gases, which might be useful in 
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investigating magnetism and superconductivity in electronic 
systems using dipolar gases as an experimental test-bed. 
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SUPPLEMENTAL MATERIALS 

1. The effective quasi-two-dimensional DDI 

The formal expression for V^p'^(r) was given in the main 
text as an integral. The integration, however, can be done 
analytically and we find V^g(r) = [D'^ / {2^/2TTal)\v{r / a^,), 
where: 

v{x) = e"!^ [(2 + x") Koix'^lA) - x^Ki{x'/A)\ , (6) 

and Kn{x) denotes the modified Bessel function of the second 
kind of order n. For small x, we find: 

vix) ^ -2\ii{x^/8) +0{x^), (7) 

while for large x, we obtain the asymptotic expansion: 

1 9 



v{x) = 2V27r 



c3 2x5 



+ 0{x-') 



(8) 



The Hartree-Fock energy density at T = can be easily 
calculated as: 



£ — £0 



1 



X 1 



4Ji(a;)- 



. 128 , 



(9) 



where Eq — fi^Trn/m in the kinetic energy of a non-interacting 
spin-polarized 2D Fermi gas, and gd and were defined in 
the main text, and I{az ) can be expressed in terms of special 
functions: 



w ^ 5 ^ /3 „ 5 5 „ o 



45 



8V2nal (Ei(87ra^) _ log(87ra2; 
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- V2 (16™^ - 3) e^™' Erf (VS^a^) + 32^2 

- 8V2 (7 - 3) nal - 24 ~ 3\/2 (e^'^"' - l) 

Expanding /(Sz) for small Gz, we get: 
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2\/2.gda, +0(5<ja2). (11) 



step, we parametrize the QMC energies of 2DDFG reported 
in Ref. 12j. Motivated by the perturbative analysis of 2DDFG 
given in Ref. O], and qualitative expectations about the de- 
pendence of the energy on g^, we parametrize the 2DDFG 
correlation energy using the following ansatz: 



.igd) = {Agj + Bgl + Cg^d) 



X log 



cigd 



3/2 

C3/25d 



C25d 



(12) 



Here, edgd) = £c(<?rf)/eoi where £0 — h^Tin/m is the kinetic 
energy of the non-interacting system. In the limit gd 1, 
the above ansatz gives £c ~ —Ag'^\og{cigd) + 0{g^J'^), 
which is in agreement with the perturbation analysis of Lu 
and Shlyapnikov |3|, who find A = —1/4 and ci « 1.43. 
Here, we do not fix these parameters using the analytic 
solution and keep both as fitting parameters. Although we 
are only concerned with the liquid phase, we remark that 
the above ansatz also assumes the correct expected form 
^ ck/V^d + /3 + 75d in the crystal phase |2|. Optimizing the 
parameters in the range gd G [0, 28], we find A — —0.2511, 
B = -0.192, C -0.01093, ci = 1.492, C3/2 = 0.004919, 
C2 = 0.5814 and C3 = 0.02382. Note that the optimized 
values of A and ci are very close to the known analytical 
values. The fit is clearly excellent. The RMSE is 5 x 10^^. 

In order to extract the kinetic energy density of 2DDFG 
from its total energy density, we use the Virial theorem. The 
general Virial theorem due to Cottrell and Paterson |4 1 implies 
2£k - jr{T,^=i r^ • dV/dn) +eds/di = 0, where V is the 
two-body interaction operator, N is the number of particles 
and £ is a linear dimension of the system size. Strictly in 2D, 
DDI is a homogeneous function of order —3 and the Virial 
term evaluates to 3 £int- Noting that gd oc i^^ , we get 2£k + 
3£int = 2£ + gddg^e[gd]- Using £ = £k + £i„t, we finally get: 



gddg^e{gd) ^ £o[l - gddg^Scigd)] 



(13) 



The pair-distribution function of 2DDFG is reported in 
Ref [2 1 for several values of gd, using which we constructed 
an analytical representation for it in the spirit of the work by 
Gori-Giorgi et al. |5|. We will report the technical details of 
this procedure in a separate work. 



3. qDFG in isotropic traps 



2. Remarks on the appHcation of the variational mapping 
method to 2DDFG 

The application of the variational mapping method to 
2DDFG is similar to mapping to 2DEG, as described in 
Ref H]. As mentioned in the main text, the only required 
information about the references system are (1) the kinetic en- 
ergy density, and (2) the pair-distribution function. As a first 



3.1. The equilibrium density profile 

The equilibrium density of the trapped gas can be obtained 
by solving Eq. (6) in the main text. For a non-interacting gas, 
£[n] £0 = Tili^n/m and we find: 

ne,,o(r) = ^(i?^P.o-r^), (14) 
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FIG. 4. The equilibrium radius (a) and density at the center of 
the trap (b) of qDFG at T = as a function of gd for various 
quasi-two-dimensionality factors (from bottom to top, \/27raz = 
0,0.1,0.2,0.3,0.4). The plots are shifted up in 0.2 increments 
for better visibility. The solid and dashed lines correspond to the 
correlated and mean-field results respectively. Note that 7?tf,o = 
{8Ny^^[h/{mujo)]^ andncq,o(0) = ma;o(2iV)^''V(27rfi^)- 



given by: 

= 2i^i(-n, n+1,1., rVi?^F,o), (15) 

corresponding to the eigenvalue Vn — ■\/2n(n + 1), where 
ri € N. In the presence of interactions, the eigenfunctions can 
no longer be found analytically. We find the spectrum of Q 
by calculating the matrix elements of Q in the non-interacting 
basis and diagonalizing it. The spectrum consists of a sin- 
gle unphysical solution, 'J ^ B^^ with eigenvalue 0. The 
oscillation frequency of the lowest lying physical monopole 
mode, rimon, corresponds to the smallest non-zero eigenvalue. 
We found that truncating the basis set to the first 10 non- 
interacting eigenfunctions yields llmon to a relative accuracy 
level of 10"^. 

3.3. Evaluating moments of the monopole response function 



where i?TF,o = (27V)i ^2 [/i/(™^o)]^/^ is the Thomas- 
Fermi radius of the cloud and N is the number of particles. 
In the presence of interactions, the differential equation for 
Ueq may only be solved numerically. We arrive at a solution 
for rieq by first finding upper and lower bounds to the density 
at center of the trap, no using a direct search, and then solving 
for the function noq(r) that satisfies the global particle number 
constraint, /d^rncq(?') = N, by bisecting the bounding in- 
terval. Convergence within a relative error tolerance of 10^^ 
is achieved usually within 20 bisections. 

Fig. [4] shows (a) the equilibrium radius of the cloud (b) the 
density at the center of the trap. The plots shows both the re- 
sults based on Hartree-Fock LDA energy density functional 
[Eq. [9) (red dashed lines) and the correlated energy obtained 
by variationally mapping to 2DDFG (black lines). Correla- 
tions make the gas more compressible and results in (1) larger 
density of particles at the center of the trap, and (2) smaller 
equilibrium radius. 



3.2. Numerical solution of the linearized hydrodynamic equation 

The linearized hydrodynamic equation for the density vari- 
ations (Eq. 7 in the main text) can be solved as follows. First, 
we note that the equation can be put in a self-adjoint form by 

changing variables to 5* = [rn/toq'T-cq] '^'^ We then 
obtain {df + G)^ ^ 0, where G = -B^^d^ ■ {n^c^^r[B^^^} 
and Scq = (mKcqnoq)^^^^- is easily verified that Q 
is self-adjoint with respect to the inner-product (^i|5'2) — 
J^^ d^r vl'*(r)*2(r), where fi+ = {r : noq(r) > 0}. More- 
over, ^ is a Hermitian operator on the functions that remain 
finite on the boundary of 51+. 

In the absence of interactions, we can find the eigenfunc- 
tions of G analytically. Here, we are interested in monopole 
oscillations and hence, only the zero angular momentum so- 
lutions are needed. The (unnormalized) eigenfunctions are 



We obtained the upper bound for flmon using the exact 
relation between the moments of the response function, i.e. 

2 (3) (1) 

(/ifimon) < "^Tw /'^TW - monopole excitation operator 
and its corresponding response functions are Sili 
and XM{i^) = h-^ J dte'-'^*{[M{t),M{0))cq respeaiwely. 
The first and third moments can be calculated using sum 
rules Q: 



(1) 
rriM = 

(3) 
n^M = 



[[M,H],M]\0)/2, 

[[[M,n],n],[n,Mm/2, 



(16) 



where |0) is the equilibrium state of the trapped gas. A direct 
calculation using commutation relations yield: 



(3) 



m 



(1) 

M 



2wn 



4(fK + ftrap) + £■ 



8£, 



trap 



(17) 



where: 



£k = y"d^rnoq(r) £KKq(r)], 

Strap = y"d^rncq(r) t/trap(r), 

= y d^r"-oq(r)£v2Kq(r)] 



(18) 



Here, [n] is given by: 



[n] 



h nn 1 



9d[nl 



(19) 



In the above equation, v'^'^\x) — xdxv(x) + d^v{x), and 
gd[n] and az[n] are the local DDI strength and quasi-two- 



dimensionality. In Eq. (19 1, the pair-distribution function is 
that of the 2DDFG that variationally maps to {gd [n],dz [n] ) . 
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